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1  INTRODUCTION  AND  OVERVIEW 


This  report  describes  an  effort  to  provide  signal  processing  tools  and  perform 
experimental  design  for  unconventional  imaging  and  parameter  estimation  from  laser- 
illuminated  targets  at  long  ranges,  using  arrays  of  detectors  without  the  need  for 
imaging  optics.  This  allows  for  an  inexpensive,  lightweight,  rapidly  steerable  receiver 
capable  of  very  fine  resolution.  This  approach  could  result  in  a  system  for  rapidly 
gathering  information  about  missiles  or  aircraft  in  flight. 

As  illustrated  in  Figure  1-1,  the  approach  is  to  illuminate  a  distant  target  with  a 
coherent  laser.  Instead  of  forming  an  image  with  bulky  imaging  optics  (a  telescope),  we 
instead  measure  the  laser  light  backscattered  from  the  target  with  an  array  of  light- 
bucket  detectors  with  no  intervening  optics.  This  measurement  scheme  allows  for  a 
much  cheaper,  lighter-weight  receiver  than  a  conventional  telescope.  The  array  could 
either  be  filled  with  detectors  or  could  be  sparsely  populated  with  detectors,  making  it 
even  more  light-weight  and  inexpensive.  The  array  could  either  be  a  heterodyne  sensor 
(mixing  the  light  from  the  target  with  light  direct  from  the  laser)  or  direct  detection 
(measuring  only  the  intensity  of  the  light  from  the  target).  In  either  case,  significant 
computer  processing  is  required  to  convert  the  output  from  the  detectors  into  usable 
information,  whether  in  the  form  of  an  image  or  in  the  form  of  useful  target 


Laser 

Figure  1-1.  Laser  Illumination  and  Detection  for  Imaging  and 
Parameter  Estimation. 
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parameters,  such  as  shape,  length,  width,  and  trajectory. 

At  the  beginning  of  the  program,  it  was  thought  that  we  might  receive  heterodyne 
data  from  another  effort,  and  that  we  could  use  it  for  imaging  experiments.  As  the 
program  developed,  it  became  apparent  that  heterodyne  data  would  not  be  available  in  a 
timely  fashion,  so  most  of  our  effort  was  aimed  a  using  direct-detection  data.  This  data, 
a  coherent  laser  speckle  pattern,  is  equivalent  to  the  intensity  of  the  Fourier  transform 
of  the  coherently  illuminated  object.  It  is  useful  primarily  for  target  parameter 
estimation. 

Under  this  effort  we  accomplished  the  following. 

We  demonstrated  the  utility  of  a  sidelobe-reduction  algorithm  and  a  superresolution 
algorithm  to  coherent  optical  imagery,  such  as  heterodyne  array  imaging.  Although  the 
algorithms  were  developed  for  synthetic  aperture  radar  imagery,  we  showed  that  they 
also  worked  well  for  an  optical  image  that  contains  fine  detail  and  specular  reflections, 
as  would  often  be  found  for  a  laser-illuminated  object.  It  gave  less  improvement  for  a 
bland,  diffuse  object.  This  is  described  is  Appendix  A. 

We  catalogued  the  target  parameters  that  could  be  characterized  by  measuring  laser 
speckle  parameters.  This  described  in  Section  2. 

We  developed  an  approach  to  achieving  improved  accuracy  in  estimating  the  shape 
of  the  object  from  laser  speckle  patterns,  using  the  sidelobe-reduction  algorithm  and  an 
autocorrelation  unfolding  approach.  This  is  described  in  Section  3. 

We  analyzed  the  signal-to-noise  ratio  for  laser  imaging  with  sparse  arrays  of 
detectors,  and  showed  the  trade-off  between  the  number  of  detectors  and  the  degree  of 
array  sparsity.  This  is  described  in  Section  4. 

We  performed  calculations  of  isoplanatic  angles  that  would  impact  any  experiments 
at  ISTEF,  as  shown  in  Section  5. 
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We  developed  a  plan  for  experiments  to  be  performed  at  ISTEF  for  measuring 
target  shapes  and  motions  from  backscattered  laser  speckle  patterns.  This  is  described 
in  Section  6. 

Section  7  gives  additional  detail  on  tracking  via  speckles,  one  of  the  target 
parameters  discussed  in  Section  2. 

Imaging  correlography  involves  averaging  the  autocovariance  of  the  pupil-plane 
laser  speckle  intensity  pattern  to  arrive  at  an  estimate  of  the  magnitude  of  the  Fourier 
transform  of  the  object’s  incoherent  reflectivity  function,  then  performing  phase 
retrieval  on  that  data  to  reconstruct  an  image.  We  developed  a  method  for 
preconditioning  autocorrelation  estimates  for  imaging  correlography.  Low-frequency 
lags  are  used  to  control  the  smoothness  of  the  reconstructed  image.  This  method  was 
demonstrated  on  laser  speckle  data  gathered  in  the  laboratory.  These  results  were 
produced  at  Michigan  Technological  University  and  are  presented  in  Appendix  B. 

A  portion  of  this  effort  is  described  in  the  presentation  and  proceedings  paper  in 

[1.1] ,  which  is  included  as  Appendix  A.  Reference  [1.2]  describes  some  of  the  results  of 
a  previous  effort  that  were  presented  in  a  conference  during  the  current  effort. 

[1.1]  “Application  of  Superresolution  Algorithm  to  Optical  Coherent  Imaging,”  J.R. 
Fienup,  M.R.  Kosek,  and  H.C.  Stankwitz,  in  Proc.  SPIE  2827-08  Digital  Image 
Recovery  and  Synthesis  III  (August  1996)  (invited  paper). 

[1.2]  “Laser  Imaging  with  Sparse  Pupil-Plane  Detectors,”  J.R.  Fienup,  to  the  Annual 
Meeting  of  the  O.S.A.,  Portland,  Oregon,  11-15  September  1995. 

References  are  at  the  end  of  each  section. 
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2  TARGET  PARAMETER  ESTIMATION 

If  we  illuminate  a  target  with  a  laser  and  detect  the  backscattered,  speckled  intensity 
pattern,  as  illustrated  in  Figure  1.1,  what  can  we  tell  about  the  target? 

If  the  measured  intensity  is  I(x,  t),  where  x  is  a  2-D  spatial  coordinate  and  t  is  time, 
then  a  valuable  quantity  to  assess  is  the  correlation  function 


ljx|2  = 


L(X  +  Ax, 
<I(X,  t)>2 


where 


I  M 

<I(x,  t)  I(x  +  Ax,  t  +  At))  =  ^  Xl(x,  tm)  I(x  +  Ax,  tm  +  At)  .  (2.2) 

Speckles  with  perfect  correlation  (spatial  and  temporal)  have  l|i|  =  1 .  When  the  speckle 
decorrelates  or  “boils,”  Ijxl  <  1. 

If  the  target  rotates  with  angular  rate  Q.,  then  the  peak  of  Ipl  translates  from  the 
origin  by  a  distance  Axp  in  the  time  interval  Atj,  given  by 

AXjy 

Q  =  ^  .  (2.3) 

We  previously  demonstrated  this  in  the  laboratory.  We  can  obtain  two  speckle  patterns 
separated  by  time  At^  either  by  (a)  reading  out  the  speckle  pattern  quickly,  using  a  fast 
detector,  or  (b)  double-pulsing  the  laser,  and  detect  the  sum  of  the  two  pulses  on  a  CCD 
camera.  A  problem  with  this  measurement  is  that  the  effective  speckle  velocity  at  the 
detector, 


Vs  =  ^  =  2QR,  (2.4) 

can  also  be  caused  by  a  transverse  translation  of  the  target  with  a  velocity  Vs/2.  The 
total  speckle  velocity  is  given  by 


Vc  =  2f2R  +  2v, 


target 
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Take  the  case  where  a  target  is  at  a  range  of  R  =  100  km  and  Q  =  0.2  radians/sec.  Then 
2£2R  =  40,000  m/sec,  which  is  much  faster  than  any  missile.  Hence  if  the  target  is 
rotating,  then  its  rotation  rate  can  easily  dominate  its  velocity  in  its  effect  on  the 
velocity  of  the  speckles.  By  observing  the  speckle  from  the  target  over  time,  it  might  be 
possible  to  determine  its  rough  trajectory  and  its  velocity  and  separate  the  effects  of  its 
rotation  rate.  Since  the  speckle  velocity  is  so  sensitive  to  the  rotation  rate,  this  affords  a 
potentially  very  sensitive  method  of  determining  the  rotation  rate  of  the  target.  Section 
7  gives  additional  details  on  measuring  rotation  rates. 

The  half-width  of  l)xl  is  related  to  the  with,  do,  of  the  object  through  the  relationship 

A.R 

^^width  =  •  (2.6) 

Uq 

Hence  we  can  determine  the  length  and  width  (projected  onto  a  plane  perpendicular  to 
the  line  of  sight)  of  the  target.  As  described  in  Section  3,  we  can  also  determine  more 
details  about  the  shape  of  the  object. 

If  the  object  consists  of  two  parts  (for  example,  a  bus  and  an  ejected  warhead),  with 
separation  S  between  them,  then  l)il  will  have  a  main  lobe  and  a  side  lobe,  with  the 
separation  between  them  being  XR/S.  Hence  we  can  measure  the  distances  between 
separated  parts  of  the  object  (projected  onto  a  plane  perpendicular  to  the  line  of  sight). 

If  the  object  is  changing,  say  due  to  thermal  effects,  then  the  speckle  pattern  will 
decorrelate  with  time,  when  the  object  is  stationary.  Hence  the  decay  of  l)j,(t)l  with  time 
is  a  measure  of  changes  in  the  object. 

If  we  measure  at  multiple  wavelengths,  then  the  decorrelation  of  l|il  with  wavelength 
gives  a  measure  of  the  depth  of  the  object  (along  the  line  of  sight). 

Several  other  discriminants  that  could  be  considered  include  the  following.  The 
color,  or  spectral  reflectivity,  of  the  object  can  be  used  to  determine  its  surface 
materials.  Polarized  light  can  be  used  to  determine  something  about  the  surface  texture 
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of  the  object  (rough  or  smooth),  and  whether  there  is  significant  3-D  geometry  effects 
(multiple  bounces  before  the  light  returns). 
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3  TARGET  SUPPORT  ESTIMATION  FROM  THE  AUTOCORRELATION 


By  the  “support”  of  the  object  we  mean  the  set  of  points  over  which  it  may  be  nonzero, 
or  its  shape  or  outline.  We  can  determine  bounds  on  the  support  of  the  object  from  the 
direct-detected  laser  speckle  pattern,  even  if  the  atmospheric  turbulence  causes  a  phase 
error  in  the  optical  fields  that  would  blur  the  image  gathered  by  a  telescope.  This  will 
allow  us  to  determine  the  length,  width,  and  orientation  of  the  object,  which  will  be 
useful  in  classification  and  trajectory  estimation.  This  section  describes  advanced 
methods  we  developed  for  support  estimating  from  laser  speckle  patterns. 

Suppose  that  the  ideal  optical  field  at  the  detector  array  is  Fk(u),  and  the 
corresponding  diffraction-limited  image  is  f^fx),  where  k  is  the  number  of  the  speckle 
realization  and  u  and  x  are  2-D  spatial  coordinates.  If  we  could  make  heterodyne 
measurements  of  the  optical  fields,  we  would  have  Gdk(u)  =  Fk(u)  exp[i<t)(u)],  where 
(t)(u)  are  the  phase  errors.  By  inverse  Fourier  transforming  we  would  arrive  at  a 
blurred  image  gdk(x).  If  instead  we  measure  the  intensity  pattern  IGdk(u)l^,  we  can 
inverse  Fourier  transform  to  arrive  at  the  autocorrelation  function 

rk(x)  =  .9‘“^[IGdk(u)|2]  =  X  gdk(x’)  gdk(x'  -  x)  =  X  fk(x')  fk(x’  -  x)  .(3.1) 

x'  x' 

Note  that  the  autocorrelation  of  the  blurred  image,  gdk(x),  is  equal  to  the 
autocorrelation  of  the  ideal  image,  fk(x),  and  it  can  be  computed  from  the  direct- 
detected  speckle  intensity  without  the  need  for  heterodyne  detection.  By  averaging  over 
many  speckle  realizations,  we  can  arrive  at  a  speckle-averaged  autocorrelation, 

ri(x)  =  K-15;  |rk(x)|2  .  (3.2) 

k 

From  this  autocorrelation  we  can  estimate  the  shape  of  the  object  using  techniques 
described  in  Reference  3.1.  Figure  3-1  illustrates  this  approach. 
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Autocorrelation  Support 


Figure  3-1.  Triple  Intersection  Approach  to  Estimate  Object  Support.  If  a  line  can  be 
drawn  tangent  to  the  support  of  the  autocorrelation  so  that  it  intersects  it  at  exactly  two 
points  (ai  and  a2),  then  the  intersection  of  three  translates  of  the  autocorrelation 
support  (centered  at  the  origin,  at  a^  ,and  at  a2)  define  a  “locator  set”  that  must  contain 
any  object  that  could  have  given  rise  to  the  autocorrelation  support  [3.1].  In  this  case, 
the  locator  set  is  exactly  the  support  of  the  original  object.  More  often  the  locator  set  is 
larger  than  the  true  object  and  gives  an  upper  bound  on  the  support  of  the  object. 


We  performed  support-estimation  experiments  with  images  of  missiles,  which  are 
germane  to  the  mission  of  SDIO.  Figure  3-2  shows  a  triple  intersection  result  assuming 
we  had  available  the  exact  autocorrelation  support  for  a  missile.  In  this  case  the  triple 
intersection  tightly  bounds  the  shape  of  the  missile,  giving  us  good  information  about 
the  length,  width,  and  shape  of  the  missile.  However,  such  idealized  autocorrelation 
information  is  not  available.  The  problem  is  that  autocorrelations  computed  from  laser- 
illuminated  objects  suffer  from  both  diffraction  and  speckle.  Diffraction  causes 
sidelobes  that  are  energy  outside  the  true  support  of  the  autocorrelation,  making  the 


8 


ERIM 

International,  Inc. 


autocorrelation  appear  larger  than  it  really  is.  Speckle  causes  energy  drop-outs  within 
the  support  of  the  autocorrelation,  making  parts  of  it  to  disappear. 


Support 

Autocorrelation  Support 


Triple  Intersection  of  Autocorrelation  Supports 

Figure  3-2.  Triple  Intersection  Approach  to  Estimate  Missile  Support. 


We  attacked  the  problems  of  speckle  drop-out  and  diffraction  to  get  better  estimates 
of  the  autocorrelation  support,  which  would  allow  us  to  get  better  estimates  of  the 
shape  of  the  object.  Figure  3-3  shows  a  block  diagram  of  the  steps  we  investigated  for 
accomplishing  this.  Figure  3-4  shows  an  example  of  each  step  for  a  simulated  image  of 
a  missile.  The  parts  (a,  b,  etc.)  of  Figure  3-4  zure  labeled  in  the  diagram  in  Figure  3-3 
and  are  referred  to  in  the  discussion  below. 

We  started  with  a  CAD/CAM  model  of  a  missile  and  computed  a  rendering  of  it  (a). 
From  that  we  computed  speckle  realizations  and  Fourier  transformed  them  to  get  the 
coherent  optical  fields  at  the  detector,  F]f(u).  The  detector  measures  only  the  intensity. 
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IPk(u)l^  (b).  We  want  to  compute  the  Fourier  transform  of  this  to  arrive  at  the 
autocorrelation  function.  However,  we  first  subtract  the  mean  intensity  to  get  rid  of  a 
delta-function  term  that  would  otherwise  appear  superimposed  on  the  autocorrelation 
function  [3.2].  To  estimate  this  mean,  we  average  the  intensity  over  the  area  of  the 
detector  (u)  or,  in  2-D  (u,  v),  and  over  the  speckle  realizations.  We  then  inverse 
Fourier  transform  the  result  to  arrive  at  a  speckled  autocorrelation  function  (c).  To 
reduce  the  sidelobes  of  the  autocorrelation  function,  we  apply  the  Spatially  Variant 
Apodization  (SVA)  algorithm  (d)  [3.3].  Then  we  average  the  intensities  of  the  speckled 
autocorrelations  over  all  the  available  speckle  realizations  to  arrive  at  a  speckle-reduced 
autocorrelation  (e,  shown  for  averaging  5  realizations).  In  the  example,  note  that  some 
small  residual  sidelobes  are  present,  emanating  from  the  origin.  This  is  due  to  an 
imperfect  estimate  of  the  bias  that  was  subtracted.  We  threshold  this  result  (e)  to  arrive 
at  an  estimate  of  the  support  of  the  autocorrelation  function  (f).  We  presently  pick  the 
threshold  by  trial  and  error.  Further  research  would  be  needed  to  pick  the  threshold 
value  automatically.  We  fill  in  the  speckle  nulls  in  this  autocorrelation  support  by  a 
“closing”  operation  from  image  algebra  (a  dilation  followed  by  an  erosion),  to  get  our 
final  autocorrelation  support  estimate  (g).  We  then  use  the  triple  intersection  rule 
described  above  to  get  an  upper  bound  on  the  object  support  (h).  We  use  this  upper 
bound  as  an  estimate  of  the  object  support.  As  suggested  in  Figure  3-4  (i),  we  can  use 
this  to  estimate  the  length,  width,  and  angle  of  the  target.  This  information  should  be 
useful  for  missile  typing  (classification)  and  determination  of  missile  trajectory.  Recall 
that  this  information  arises  from  the  pupil-plane  intensity  patterns,  as  shown  in  Figure 
3-4(b),  which  looks  nothing  like  an  image  and  does  not  require  optics  to  form  an 
image. 
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Multiple  aperture  plane 
speckle  realizations 
Fk(u,v) 


(c)  (d) 


Object 

/px  autocorrelation 
'  ^  estimate 


To  shape-based 
analysis  and 
identification 


Figure  3-3.  Block  Diagram  of  Steps  to  Improve  Image  Support  Estimation. 
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(a)  Incoherent  image 


(d)  AC  after  SVA  (e)  Average  AC 


(c) 


Coherent  Autocorrelation 


(g)  Filled  AC  Support  (h)  Object  Support  Bound  (i)  Potential  Processing  Result 
Figure  3-4.  Image,  Fourier,  and  Autocorrelation  Steps  in  Estimating  Object  Support. 
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Figure  3-5  shows  the  results  of  averaging  different  numbers  of  speckle  frames, 
from  K  =  1  frame  (no  averaging),  to  K  =  <»,  which  is  equivalent  to  an  incoherent 
autocorrelation.  The  object  support  estimate  improves  as  we  increase  the  number  of 
speckle  frames  averaged. 

With  these  results  we  have  shown  that,  with  advanced  signal  processing  methods,  we 
can  obtain  reasonable  estimates  of  the  length,  width,  and  orientation  of  a  missile  from 
the  intensity  of  laser  light  backscattered  from  it,  without  using  imaging  optics.  These 
simulations  assumed  that  the  laser  power  was  high  enough  that  speckle  noise  statistics 
dominates  over  photon  noise.  Further  research  would  be  necessary  to  determine  the 
effect  of  photon  noise  on  the  result. 
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Figure  3-5.  Autocorrelation  and  Object  Support  Estimates  for  Different  Numbers  of 
Speckle  Realizations. 
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4  ARRAY  SPARSITY  AND  SPECKLE  REALIZATIONS 


The  signal-to-noise  ratio  (SNR)  of  our  autocorrelation  estimates  computed  from  the 
pupil-plane  intensity  measurements  from  a  laser-illuminated  object  depends  on  the  light 
level,  the  size  of  the  detector  array,  the  sparsity  of  the  array  (the  array  fill  factor),  and 
the  number  of  speckle  realizations  over  which  we  average.  Suppose  that  we  pick  an 
array  diameter  to  achieve  the  desired  resolution  on  the  target  (p  =  X,R/D,  where  p  is 
resolution,  X  wavelength,  R  range,  and  D  diameter  of  the  detector  array).  We  can  make 
the  array  cheaper  and  more  light  weight  by  sparsely  populating  it  with  detector 
elements,  but  then  we  require  a  greater  number  of  speckle  realizations  (i.e.,  a  greater 
collection  time)  to  achieve  the  required  SNR. 

For  heterodyne  array  imaging,  the  SNR  at  a  given  spatial  frequency  (u,  v)  is  given 
by  [4.1] 

SNR(u,  v)  «  K  ngOTF(u,  v)  lp,(u,  v)l  ,  (4.1) 

and  for  correlography  (measuring  speckle  intensities  in  the  pupil  plane)  [4.2], 

SNR(u,  V)  -  ^(K/3)  ngOTF(u,  v)  lp(u,  v)l^  ,  (4.2) 


where 


K 

“s 

ngOTF(u,  V) 

l|i(u,  v)P 


=  number  of  speckle  realizations 
=  number  of  independent  pupil-plane  detectors 
=  number  of  redundant  measurements  at  (u,  v) 

=  normalized  power  spectral  density  of  object  intensity 


The  approximations  above  are  good  when  speckle  noise  dominates  (for  more  than  two 

detected  photons  per  speckle)  and  for  the  mid-  to  upper  spatial  frequencies  where  lp,(u, 
2 

v)l  is  small. 


We  see  that,  for  either  modality,  the  SNR  goes  as  the  square  root  of  ngOTF(u,  v), 
which  in  turn  depends  on  the  number  of  detectors  in  the  array,  the  sparsity  of  the 
array,  and  the  particular  configuration  (geometrical  arrangement)  of  the  elements  of 
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the  array.  A  decrease  in  ngOTF(u,  v)  by  some  factor  can  be  made  up  for  by  a 
proportional  increase  in  K. 


2 

The  power  spectral  density,  l|i.(u,  v)l  ,  of  the  object  depends  greatly  on  the  type  of 
object,  as  illustrated  in  Figure  4-1.  For  point-like  objects  (e.g.,  stars)  it  falls  off  slowly 
with  spatial  frequencies.  As  the  object  becomes  more  extended,  l|j,(u,  v)l  drops  off 
more  rapidly  with  spatial  frequency.  A  rule  of  thumb  (by  F.  Roddier)  is  that  if,  at  a 
given  maximum  spatial  frequency  %  we  can  resolve  n^p  bright  points,  then  l|J,(Uj.)l  « 
l/n^,p(Uj.).  Hence  the  performance  depends  on  the  spatial  frequency  content,  or  the 
stmcture,  of  the  object,  as  well  as  on  the  amount  of  signal  collected. 


Figure  4-1.  Object  Power  Spectral  Densities 
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2 

We  find  that,  since  l|j,(u,  v)l  falls  off  for  the  highest  spatial  frequencies  as  does 
OTF(u,  v),  we  can  rarely  have  good  SNR  at  the  highest  spatial  frequencies.  We  are 
better  off  comparing  performance  at  the  middle  spatial  frequencies.  Hence  we 
concentrate  on  the  spatial  frequency  u  =  (1/2)  Uj^^x- 

Now  consider  the  term  n30TF(u,  v),  the  number  of  redundant  measurements.  As 
illustrated  in  the  top  half  of  Figure  4-2,  the  OTF,  which  is  given  by  the  autocorrelation 
of  the  aperture  function,  normalized  to  be  unity  at  the  origin,  is  about  0.5  for  a  filled 
aperture  at  u  =  (l/2)Ujjjax-  That  is,  when  the  aperture  is  sheared  by  half  its  width,  about 
half  of  the  sheared  elements  fall  on  the  positions  of  the  original  elements,  and 
OTF(Uinax/2)  =  1/2.  For  a  nonredundant  aperture  (not  shown,  since  none  are  known  to 
exist  for  the  2-D  aperture  sizes  shown),  there  would  be  only  a  single  overlapping 
element,  and  so  the  OTF  would  be  1/N,  where  N  is  the  number  of  nonzero  elements 
and  N/M^  is  the  aperture  fill  ratio.  So  nonredundant  apertures  have  very  low  values  of 
n50TF(u,  v)  and  require  M^/N  times  more  speckle  realizations  that  a  filled  aperture  of 
the  same  diameter.  The  bottom  example  in  Figure  4-1  shows  an  aperture  with  N  =  54 
subapertures  designed  by  P.  Henshaw  [4.3]  to  be  minimally  redundant.  Depending  on 
the  spatial  frequency,  it  has  n50TF(u,  v)  =  anywhere  between  one  and  several  overlaps, 
and  has  6  for  the  spatial  frequency  shown,  where  OTF(u,  v)  =  6/54  «  0.1 1.  Figure  4-3 
shows  the  arrow  and  box+2  array  patterns,  which  we  invented.  The  former  has  ngOTF 
>  1  and  the  latter  has  ngOTF  >  2  for  all  spatial  frequencies,  using  a  relatively  small 
number  of  array  elements,  and  both  arrays  are  easily  scalable  to  any  size  desired. 
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Genetic  Algorithm  --  54  subapertures 

Figure  4-2.  Array  Patterns  and  their  OTF's.  Top:  filled  aperture;  bottom:  pattern 

generated  by  genetic  algorithm 
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“Arrow”  Array 
3M  -  3  subapertures 
to  cover  same  area 
as  M  X  M  filled  array 
with  no  Fourier  gaps 


“Box  +  2”  Array 
4M  -  2  subapertures 
to  cover  same  area 
as  M  X  M  filled  array 
with  minimum  of 
double  redundancy 


Figure  4-3.  Arrow  and  Box+2  Array  Patterns. 


Now  let’s  compare  a  sparse  array  with  N  elements  to  a  filled  array  with  MxM  = 
elements.  Both  cover  a  (2M  -  1)  x  (2M  -  1)  array  in  Fourier  (OTF)  space.  Including 
the  redundancy  factor  for  each  spatial  frequency,  with  N  elements  we  measure  2N(N- 
1)  frequencies  (not  counting  d.c.).  An  ideal  (fictitious)  nonredundant  array  covering 
(2M  -1)  X  (2M  -1)  Fourier  samples  would  require  N  elements,  where  2N(N  -  1)  = 
(2M  -  1)2  -1,  or  N  «  V2  M  for  large  M.  A  realizable  minimally  redundant  array,  the 
arrow  array,  requires  N  =  3M  -  3  elements.  The  cost  savings  in  array  elements, 
compared  with  the  filled  aperture,  is  a  factor  of  M2/(3M  -  3)  «  M/3,  which  can  be  a 
large  number.  But  at  u  =  (l/2)Ujjjax»  the  filled  aperture  has  UgOTF  =  M2/2,  versus  1  for 
the  weakest  spatial  frequency  of  the  minimally  redundant  array.  Hence  the  factor  of 
M/3  in  cost  savings  must  be  weighed  against  the  factor  of  M2/2  in  integration  time  to 
get  that  factor  times  as  many  speckle  realizations. 
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The  number  of  speckle  realizations  that  are  available  limits  how  great  an  increase  in 
SNR  is  possible.  To  get  a  new  speckle  realization,  the  target  must  rotate  through  an 
angle  D/(2R)  =  X/(2p),  where 

D  =  aperture  diameter, 

R  =  range  to  target, 

A,  =  wavelength,  and 
p  =  resolution. 

However,  if  an  object  rotates  too  much,  it  no  longer  looks  like  the  same  objects.  We  use 
a  rule  of  thumb  that  the  object  can  rotate  by  p/d^  radians  before  it  changes 
significantly,  where  d^  is  the  diameter  of  the  object.  This  allows  for  (p/do)(2p/A.)  = 
2p/(doA,)  speckle  realizations.  For  example,  if  the  object  diameter  is  d^  =  5  m,  p  =  0.20 
m,  and  A,  =  1  pm,  then  we  can  have  up  to  16,000  speckle  realizations,  if  the  laser  and 
detector  are  fast  enough  to  detect  in  the  time  that  it  takes  the  speckle  pattern  to  move  by 
D. 


To  achieve  a  given  SNR,  we  find  that  for  heterodyne  array  imaging  the  number  of 
speckle  realizations  we  need  is 


K 


SNR2 

ngOTF  lp|2 


(4.3) 


and  for  correlography  is 


_  3  SNR2 

“  ngOTF  Ipl4  • 


(4.4) 


If  we  assume  that  we  need  a  SNR  =  3  at  the  spatial  frequency  u  =  (1/2)  Uj^ax^  then  Table 
4.1  compares  the  number  of  speckle  realizations  for  (a)  filled  versus  nonredundant  (or 
minimally  redundant)  arrays,  (b)  for  heterodyne  imaging  and  correlography,  (c)  for 
three  different  object  types,  and  (d)  for  array  sizes  8x8,  32x32,  and  128x128.  From 
this  we  see  that  nonredundant  and  minimally  redundant  arrays  are  generally  suitable 
(require  less  than  10^  speckle  realizations)  for  coherent  illumination  only  for  point-like 
object  or  small  extended  objects.  For  large,  extended  objects,  nonredundant  arrays  are 
not  suitable,  and  one  needs  arrays  with  large  amounts  of  redundancy  (such  as  a  filled 
array). 
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Table  4.1.  Number  of  Speckle  Realizations  Required: 


object  type: 

two  glints 

small  extended 

large  extended 

MxM 

Modality 

lUl  =  0.5 

II 

p 

filled  nonred. 

filled  nonred. 

filled  <-»■  nonred. 

8x8 

heterodyne 

1^36 

28  ^  900 

n/a 

8x8 

correlography 

14  ^  432 

8,438  ^  2.7x10^ 

n/a 

32x32 

heterodyne 

1  ^36 

2  <->-900 

176  ^  90,000 

32x32 

correlography 

1  ^432 

527  ^  2.7x10^ 

5.3x10^  ^  2.7x10^ 

128x128 

heterodyne 

1  ^36 

1  ^900 

11  ^  90,000 

128x128 

correlography 

1  ^432 

33  ^  2.7x10^ 

3.3x10^  ^  2.7x10^ 
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5 ISOPLANATIC  CONSIDERATIONS 


Introduction 

We  are  interested  in  determining  what  types  of  information  can  obtained  with  a 
detector  array  acquisition  system  that  captures  a  speckle  pattern  formed  from  a  missile 
illuminated  with  coherent  laser  light.  Under  consideration  is  using  the  speckle  for  imag¬ 
ing,  target  shape  discrimination  (missile  typing)  via  speckle  shapes,  and  target  rotation 
rate  determination  via  speckle  motion.  Atmospheric  turbulence  can  affect  the  detected 
speckle  pattern  and  hence  the  quality  of  the  target  measurements.  The  most  stringent 
requirements  on  the  atmospheric  turbulence  for  the  various  techniques  is  for  generating  an 
image.  It  requires  that  the  received  wavefront  from  each  scatterer  on  the  target  propagate 
through  nearly  the  same  atmospheric  distribution.  This  is  known  as  the  isoplanatic 
requirement.  Our  goal  here  is  to  calculate  some  of  the  limitations  for  missile  imaging  via 
speckle  as  given  by  isoplanatic  relationships. 

Atmospheric  turbulence  causes  spatially  and  temporally  random  fluctuations  in  the 
index  of  refraction  of  the  atmosphere,  which  then  results  in  random  spatial  and  temporal 
variations  in  the  phase  of  an  propagating  optical  wavefront.  Although  the  fundamentally 
local  change  is  phase,  further  (far  field)  propagation  of  the  phase  aberrated  wave  can  also 
lead  to  amplitude  variations  in  the  wavefront.  In  our  considerations  here,  we  are  assuming 
that  only  phase  variations  in  the  wavefronts  are  significant  and  the  wavefronts  are  uniform 
in  amplitude  for  both  the  target  illumination  wavefront  and  in  the  scattered  wavefront. 

We  also  need  to  distinguish  between  short  exposures  recordings  of  the  speckle  pat¬ 
tern  and  long  exposure  recordings.  By  short  exposure  recording  we  mean  that  the  atmo¬ 
sphere  can  be  assumed  to  be  frozen  in  its  configuration  during  the  exposure.  T3q)ically,  this 
is  true  for  exposures  <10  ms.  Long  exposure  recordings  must  take  into  account  the  tempo¬ 
ral  changes  in  the  atmosphere.  Operationally  speaking,  we  are  considering  the  case  where 
the  target  is  illuminated  with  a  long  train  of  short  light  pulses.  For  example,  we  are  consid¬ 
ering  using  the  ISTEF  Long  Pulse  Laser  that  emits  a  continuous  series  of  0.6  jxsec  pulses 
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at  a  pulse  rate  of  3  per  second.  A  recording  of  the  speckle  pattern  for  any  one  pulse  would 
correspond  to  a  short  exposure,  while  long  exposure  recording  is  obtained  with  integration 
over  two  or  more  pulses. 


Isoplanatic  Expressions  and  Calculations 

For  speckle  imaging  we  require  that  the  targets  angular  extent  not  exceed  the 

isoplanatic  angle  Qq.  Hence,  we  have  the  relationship 


L<RdQ  (5.1) 

where,  L  is  the  target’s  largest  dimension  and  R  is  the  range  to  the  target.  The  short  expo¬ 
sure  and  long  exposure  cases  have  different  isoplanatic  angles. 

The  formula  for  calculating  Qq  is  given  by  [5. 1-5.2] 

3 

0 

2 

where  is  the  structure  parameter  of  the  atmospheric  index  of  refraction  fluctuations.  It 
is  typically  a  function  of  altitude  and  climatic  conditions,  and  is  constant  in  the  horizontal 
direction  over  uniform  terrain.  At  ground  level  we  can  expect  values  of  10'^^  to  10"*^  and 

at  an  altitude  of  10  km  values  of  10'*^  to  10'^^  for  . 

2 

For  constant  C„  (i.e.,  typically  for  horizontal  propagation)  Eq.  2  gives 

8_,_3 

2  5  6  8  6 

^  =  0.105  C/  Z?  ^  (5.3) 

X^\ 

The  overall  effect  on  the  system  response  can  often  best  be  calculated  through  the 
optical  transfer  function  (OTF).  The  atmospheric  OTF  contains  the  factor 


MTFq^  =  exp 


(5.4) 
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Hence,  at  the  isoplanatic  angle  the  OTF  is  reduced  by  e  and  the  speckle  size  is  correspond¬ 
ingly  increased. 

Also  of  interest  is  what  we  will  call  the  isoplanatic  patch  diameter  Pq  =RQo,  which 

2 

can  be  written  for  constant  C„  as 

fl 


Pq  =  0.105  • 


■^2  5 

^n'  2 

A,  -I 


(5.5) 


A  key  point  to  see  is  that  the  isoplanatic  patch  diameter  decreases  with  increasing  range. 

Hence,  there  is  a  maximum  range  for  a  desired  patch  diameter. 

2 

The  C„  constant  varies  with  altitude  and  climatic  conditions.  At  ground  level  we 
can  expect  values  of  10'^^  to  lO'^"^  and  at  an  altitude  of  10  km  values  of  10"^^  to  10'^^.  Fig- 
ures  la-lc  show  plots  of  the  isoplanatic  angle  versus  range  for  different  constant  and 
different  wavelengths.  Figure  5-2a-5-2c  show  the  corresponding  isoplanatic  patch  diame¬ 
ter  as  a  function  of  range.  It  can  be  seen  that  at  ground  level  the  isoplanatic  length  can  be 
quite  limiting.  Working  at  longer  wavelengths  and  smaller  C„  (higher  altitudes)  alleviates 
the  problem.  Figure  5-3  shows  the  critical  values  of  and  range  R  using  equation  5.5 
when  Pq  is  set  to  12  meters.  The  12  meters  is  appropriate  for  the  length  of  scud  missiles. 
Note  that  if  we  were  to  work  at  a  range  of  100  km  at  a  wavelength  of  0.5  pm  it  is  seen  that 
an  isoplanatic  patch  size  of  12  m  requires  a  of  about  10'^%  which  typically  corre¬ 
sponds  to  an  altitude  of  over  25  km.  This  altitude  is  quite  high  and,  for  example,  exceeds 
the  13  km  ceiling  of  the  AWACS. 


We  also  know  from  guide  star  work  that  the  isoplanatic  angle  can  be  a  few  arcsec- 
onds  for  zenith  viewing  at  a  wavelength  of  0.5  pm.  This  constant,  range  independent 
isoplanatic  angle,  is  due  to  the  stars  being  outside  the  atmosphere  and  the  calculations  are 
for  looking  through  the  entire  thickness  of  the  atmosphere.  Note  that  5  arcseconds  equals 
about  24  prad.  This  isoplanatic  angle  corresponds  to  a  range  of  500  km  for  a  12  m  patch 
size.  Thus,  the  target  must  be  over  500  km  away  to  be  entirely  within  the  isoplanatic  angle 
for  zenith  viewing.  This  height  exceeds  the  operating  altitude  of  any  tactical  missile  and 
hence  the  technique  is  not  useful  for  up-looking  systems. 
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As  a  point  of  interest,  we  note  that  the  form  of  the  expression  for  the  atmospheric 
coherence  length  Tq  (Fried’s  coherence  radius  [5. 1,5.2])  is  very  similar  to  Eq.  5.2  and  that 
we  can  write  Tq  =  1.76  •  Pq  for  parallel  plane  illumination  and  a  constant  C^.  Thus, 
within  a  factor  of  two,  Tq  must  be  about  the  same  size  as  the  maximum  target  size. 
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Path  Length  (meters) 


Figure  5- la.  Isoplanatic  Angle  vs.  Path  Length  at  0.5  jxm  Wavelength  for  Various 
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Figure  .5- lb.  Isoplanatic  Angle  vs.  Path 
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6 ISTEF  EXPERIMENTAL  DESIGN 


The  objective  of  this  proposed  plan  is  to  experimentally  demonstrate  at  BMDO/ 
ISTEF  (Ballistic  Missile  Defense  Organization  Innovative  Science  and  Technology  Exper¬ 
imental  Facility)  the  robustness  and  sensitivity  of  speckle  based  target  shape  discrimina¬ 
tion  and  motion  determination  of  remote  targets.  We  propose  to  use  ISTEF  test  ranges, 
other  facilities,  and  available  equipment  as  much  as  possible. 

Concept 

The  basic  concept  is  to  illuminate  a  target,  such  as  a  missile,  with  a  one  pulse  or 
train  of  coherent  laser  Ught  pulses  and  to  directly  record  resulting  speckle  patterns  with  a 
detector  array.  The  conceptual  experimental  arrangement  is  illustrated  in  Figure  6-1. 


Figure  6-1.  The  Speckle  Based  Concept. 

Notice  that  no  attempt  is  made  to  optically  form  an  image  of  the  target,  but  instead  the 
intensity  (speckle)  pattern  that  is  formed  from  the  light  scattered  from  the  target  is 
recorded  with  an  electronic  camera.  By  processing  the  recorded  speckle  patterns,  target 
shape  information  and  target  motion  information  is  obtained.  The  shape  of  the  speckles 
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allows  us  to  obtain  information  about  the  shape  and  orientation  of  the  target.  The  speckle 
motion  allows  us  to  gather  information  about  the  target’s  trajectory  and/or  rotation  rate. 

The  advantage  of  using  such  simple  (non-imaging)  systems  is  that  detector  arrays 
could  be  mounted  conformally  on  the  surface  of  an  acquisition  platform  without  requiring 
expensive,  heavy,  and  bulky  optics.  Also,  the  processing  computational  load  required  to 
obtain  the  desired  information  is  small  enough  to  allow  real-time  operation. 

Selection  of  Experiments 


There  are  a  number  of  possible  speckle  experiments  listed  in  Table  6. 1  that  can  be 

Table  6. 1  Potential  Experiments 


Laser 

Detector 

Measurements 

Rotation 

Rate/ 

Direction 

2-D 

Shape 

3-D 

Shape 

Separation 

Vibration 

cw 

CCD-SF-G 

X 

X 

X 

CW 

CCD-MF-G 

X 

X 

X 

X 

CW-Tunable 

CCD-MF-G 

X 

X 

X 

Short  Pulse 

ICCD-SF 

X,E 

X,E 

Double 

Short  Pulses 

CCD-SF 

X 

X 

Long  Pulse 

CCD-SF 

X,E  ,  / 

X,E 

X,  E 

Long  Pulse 

Double  Gated 
IGCD-SF 

X,  E  ■  ■ 

Double  Long 
Pulses 

CCD-SF 

X 

Key: 

SF  -  Single  Frame,  MF  -  Multiple  Frames,  G  -  Gated 
CCD  -  Charge  Coupled  Detector  Array  Camera 
ICCD  -  Intensified  CCD 
X  - 1  km  Range,  E  -  Eastern  Range 

carried  out  based  on  various  types  of  lasers,  cameras,  and  ranges.  The  highlighted  table 
cells  are  experiments  that  can  be  accomplished  with  known  ISTEF  equipment  on  either 
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the  ISTEF  1  km  laser  test  range  (denoted  by  X  in  Table  1)  or  the  Air  Force  Eastern  Range 
(denoted  by  E  in  Table  6.1)[6.1].  The  1  km  test  range  starts  at  the  main  site  and  extends  to 
the  north.  The  Air  Force  Eastern  Range  is  a  nearby  facility  where  the  Air  Force  launches 
rockets  from  Cape  Canaveral  Air  Station  and  at  which  laser  imaging  has  been  accom¬ 
plished  at  ranges  up  to  5  to  7.5  km  [6.2]. 

There  are  basically  3  different  measurements  (rotation/motion  of  a  target,  2-D  pro¬ 
jected  shape  of  a  target,  and  separation  between  targets)  that  could  be  accomplished  at  the 
two  ranges.  Properties  of  the  two  available  ISTEF  lasers  of  interest  are  given  in  Table  6.2 
and  properties  of  two  available  cameras  of  interest  are  given  in  Table  6.3  [6.1]. 


Table  6.2  Lasers 


Parameter/Laser 

Short  Pulse 

Long  Pulse 

Type 

Nd:YAG 

Nd:  Glass  Slab 

Pulse  Length  (nsec) 

10 

600 

Pulse  Rep  Rate  (Hz) 

10 

3 

Energy  per  Pulse  (J) 

1.5  at  1064  nm 

0.5  at  532  nm 

30  at  1054  nm 

16  at  527  nm 

Table  6.3  Cameras 


Parameter/Camera 

Sony  (2)  SSC-D7 

Xybion  350  (1) 

Detector 

Silicon 

Gen  2  Intensified 

Pixel  Size  at  CCD  (|im) 

11.5  X  13.4 

1 1.5  X  13.4 

Effective  Pixel  Size  at 
Intensifier  Face  (|im) 

Not  Applicable 

17.2x20.1 

Pixel  Count 

768x493 

768x493 

Effective  Aperture  (mm^) 

8.8  X  6.6 

12.7  X  9.6 

Intensifier  Gain 

Not  Applicable 

300  to  18,000 

Minimum  Gating  Time  (ns) 

Not  Applicable 

25 

Wavelength  Range  (nm) 

0.4-1. 1 

0.48-0.88 
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Frame  Processor  & 

Grabber  Controller 


Figure  6-2.  The  Baseline  Experimental  Arrangement. 


Baseline  Experiments 

The  baseline  experimental  arrangement  is  illustrated  in  Figure  6-2.  A  full  size  mis¬ 
sile  or  surrogate  (scaled)  target  is  mounted  on  a  pedestal  such  that  it  can  be  rotated  in  vari¬ 
ous  ways  to  present  dynamic  changes  in  perspective  as  the  experimental  data  is  acquired. 
A  train  of  coherent  laser  hght  pulses  illuminate  the  target.  Scattered  Ught  from  the  target 
for  each  pulse  is  captured  as  a  speckle  pattern  by  a  digital  camera  and  acquisition  system. 
The  relative  shape  of  the  target  and  its  motion  is  obtained  by  digitally  processing  the  cap¬ 
tured  train  of  speckle  patterns. 

The  range  could  cover  short  distances  (even  indoors)  with  scaled  models  to  long 
distances  with  full  size  targets.  It  is  also  possible  that  actual  launched  missiles  could  be 
tracked  and  detected  by  modifying  current  ISTEF  laser  imagers  that  have  been  used  for 
imaging  the  NASA  Cassini  launch  [6.2].  However,  our  proposed  baseline  experimental 
setup  is  for  a  full  size  (~  0.9  m  by  10  m)  “scud  class”  of  missile  target  mounted  at  a  range 
of  1  km  using  the  ISTEF  1  km  range.  This  allows  the  speckle  based  concepts  to  be  evalu¬ 
ated  on  full  sized  targets  under  controlled  geometrical  conditions. 
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Source/Transmitter 

The  baseline  light  source  is  the  ISTEF  Long  Pulse  Laser.  It  is  a  Nd:  Glass  slab  laser 
that  puts  out  30  Joules/pulse  at  1054  nm  or  16  joule/pulse  at  527  nm  with  a  pulse  length  of 
600  ns  and  with  a  pulse  repetition  rate  of  up  to  3  pulses  per  second.  The  output  at  the  visi¬ 
ble  wavelength  (527  nm)  is  the  baseline  wavelength  of  operation.  The  pulses  are  short 
enough  (i.e.,  <  1  ms)  to  ensure  that  the  atmospheric  turbulence  can  be  considered  to  be  fro¬ 
zen  over  the  exposure  time  associated  with  each  pulse.  Each  pulse  is  to  be  spatially  and 
temporally  coherent  with  itself  such  that  the  entire  3-D  target  is  coherently  illuminated 
during  a  pulse.  However,  no  coherence  between  pulses  is  expected  since  the  pulse  rate  is 
too  slow  for  multiple  pulses  during  the  atmospheric  freeze  time. 

The  Long  Pulse  Laser  can  be  fed  into  the  Contraves  Graz  optical  tracker  mount 
installed  in  the  tower  and  dome  at  the  northeast  comer  of  the  ISTEF  site.  The  mount  has  a 
20”  reflecting  telescope  operating  at  f/2  which  can  operate  from  the  UV  through  MWIR 
wavelengths.  It  has  a  2  inch  clear  aperture  Coude  feed  which  has  been  used  to  give  532nm 
and  1064nm  illumination  for  collimated  output.  To  illuminate  the  entire  target  the  feed 
will  need  to  be  modified  to  give  an  expanding  beam.  However,  even  partial  illumination  of 
the  target  could  produce  useful  information  for  the  target  rotation  experiments. 

The  Long  Pulse  Laser  has  the  potential  to  be  double  pulsed,  which  could  give  our 
experiments  more  versatility,  but  this  double  pulsing  modification  has  not  been  imple¬ 
mented. 

Target 

We  propose  a  target  for  the  baseline  experiments  that  is  a  full-scale  mock-up  of  the 
scud  missile  type,  or  possible  an  actual  missile.  This  type  of  missile  has  a  range  of  sizes 
and  characteristics  as  listed  in  Table  6.4  [6.2].  For  analysis  and  calculation  purposes,  the 
nominal  baseline  experimental  dimensions  for  the  missile  are  taken  to  be  0.9  m  in  diame¬ 
ter  by  10  m  in  length. 
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Property  Values 

Designation  SS-1(A,  b,  B,  c, 

or  C) 

Length  9.45  to  11.58 

Diameter  0.76  to  0.914 

Length/Diameter  8 . 3 : 1  to  1 4 : 1 

Ratio 

Max  Speed  1,100  to  1,500  m/sec 

(Cut-off  Speed) 

Range  160  to  450  km 

Wingspan/Fin  0  to  2.9  m 

The  initial  setup  stage  would  be  to  mount  the  missile  in  a  vertical  position  on  a 
rotating  table  as  illustrated  in  Fig.  6.3a.The  table  should  be  elevated  to  minimize  ground 
turbulence.  The  missile  could  be  rotated  on  the  table  to  simulate  a  spinning  missile.  Parts 
of  the  missile,  e.g.  the  tails,  could  be  changed  to  different  sizes  to  see  if  the  speckle  tech¬ 
niques  can  distinguish  the  changes.  Multiple  targets  could  also  be  mounted  and  illumi¬ 
nated  simultaneously  to  see  if  they  can  be  sensed  by  the  speckle  techniques. 

The  second  stage  is  where  the  target  is  mounted  onto  a  structure  that  allows  three 
degrees  of  dynamic  rotational  freedom  in  order  to  allow  us  to  simulate  a  missile  trajectory 
except  for  the  translation.  One  rotational  axis  allows  the  missile  to  be  rotated  about  its 
length  axis  at  a  rotation  rate  of  up  to  1  revolution  per  second  with  up  to  5  degrees  of  offset 


Units 

US  Designations 

m 

m 
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to  simulate  wobble  as  might  occur  in  actual  flight.  Another  rotational  axis  (vertical)  would 
allow  the  missile  to  turn  slowly  towards/away  from  the  line  of  sight  directions  and  the 
another  (horizontal)  would  allow  it  to  slowly  be  orientated  with  different  inclinations  with 
respect  to  the  ground  plane. 

One  way  of  achieving  this  type  of  mounting  is  to  attach  the  missile  to  two  posts, 
one  at  each  end  of  the  missile.  As  illustrated  in  Fig.  6-2b  ,  the  posts  are  mounted  on  a  tum- 


(a)  Initial  Mount  (b)  More  Versatile  Mount 

Figure  6-3.  Target  Mount  Assemblies. 


table  to  allow  the  entire  missile  to  rotate  about  a  vertical  axis  Allowing  one  of  the  post  to 
be  extendable  would  give  the  inclination  rotation.  Rotating  the  entire  missile  about  its 
attachment  to  the  post  positions  would  give  main  rotation  about  its  length  axis  and  if  the 
attachment  point  to  one  of  the  posts  is  offset  from  the  missile  axis  then  a  wobble  action 
would  be  generated. 

We  note  that  a  concrete  pad  with  power  exists  at  the  end  of  the  1  km  range.  How¬ 
ever,  any  mounts  for  the  targets  will  need  to  be  constructed. 

Detector/Receiver 

The  basic  reception  process  is  to  detect  the  speckle  patterns  with  a  detector  array 
camera  and  capture  the  camera  data  with  a  digital  frame  grabber  acquisition  system.  Con- 
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ceptually,  no  lenses  are  necessary  since  the  speckle  pattern  at  the  camera  array  is  the 
desired  quantity  to  be  acquired.  However,  since  we  require  at  least  two  detectors  per 
speckle  for  each  of  the  two  orthogonal  directions,  a  lens  system  may  be  needed  to  image 
the  speckle  pattern  on  the  detector  array  with  correct  magnification  to  properly  optinuze 
sampling  the  speckle  pattern.  Further,  the  optics  can  be  used  to  restrict  the  field-of-view  to 
help  reduce  background  light  detection  and  hence  improve  the  signal-to-noise  ratio.  A 
bandpass  filter  matching  the  laser  output  wavelength  is  also  required  to  exclude  sunlight 
and  other  extraneous  light  in  the  field-of-view.  A  polarization  filter  may  also  be  added  to 
study  polarization  influences  on  the  techniques. 

For  our  baseline  experiments  we  plan  to  detect  the  speckle  pattern  using  the  ISTEF 
Xybion  350  Camera  with  its  own  mount  near  the  support  building  such  that  the  indoor 
acquisition  and  processing  equipment  can  be  used.  The  camera  is  an  intensified  CCD 
detector  type  with  768X493  pixels  and  with  an  aperture  size  of  12.7  mm  by  9.6  mm.  It  is  a 
gated  camera  that  can  be  operated  to  expose  only  during  the  time  the  laser  pulse  is  return¬ 
ing  from  the  target.  Thus,  the  background  light  during  other  times  than  during  the  expo¬ 
sure  will  not  be  detected.  Appendix  C  shows  a  MathCad  calculation  sheet  with  the  basic 
equations  and  parameters  for  the  baseline  experiments.  With  the  indicated  parameters  we 
find  that  the  number  of  photons  per  speckle  is  about  6900,  which  is  independent  of  range 
since  the  angular  speckle  size  is  only  dependent  on  the  target  size  and  operational  wave¬ 
length  and  since  we  assume  the  illumination  beam  is  adjusted  in  its  width  to  keep  the  same 
light  energy  on  the  target.  We  will  use  the  camera  directly  with  no  magnification  optics.  At 
a  1  km  range  we  have  for  the  11:1  aspect  ratio  target  about  3  detectors  per  speckle  in  one 
dimension  and  about  30  in  the  orthogonal  dimension.  This  arrangement  assumes  that  the 
camera  aperture  with  its  12.7  mm  width  (768  pixels  X-direction)  is  orientated  along  the 
short  (0.9  m)  dimension  of  the  target  and  the  9.6  mm  width  (493  pixels  Y-direction)  is  ori¬ 
entated  along  the  long  dimension  (11  m)  of  the  target.  With  this  arrangement  the  calcula¬ 
tions  in  Appendix  C  give  that  there  are  more  than  32  photons  per  detector  element  for  each 
laser  pulse.  The  camera  intensifier  can  give  a  gain  of  over  15,0(K3.  Hence,  a  large  number 
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of  photoelectron  events  and  hence  significant  output  voltage  per  pixel  can  be  generated 
even  for  a  few  laser  photons  per  pixel. 

An  additional  advantage  of  using  the  Xybion  camera  is  that  it  can  be  gated  twice 
during  the  laser  pulse  time  of  600  ns.  Thus,  if  the  target  is  rotating  we  can  obtain  two 
exposures  before  the  atmosphere  has  changed  significantly  and  demonstrate  how  the 
speckle  translation  can  be  used  to  calculate  the  motion  of  the  target.  Given  the  arrange¬ 
ment  described  above  with  no  magnification,  we  have  about  22  speckles  across  the  aper¬ 
ture  in  the  X-direction.  If  we  use  two  25  ns  exposures  500  ns  apart  and  a  maximum 
speckle  translation  of  1 1  pixels,  we  can  measure  a  rotation  rate  of  up  to  about  1  revolution 
per  second  for  the  target.  The  smear  is  only  about  half  a  pixel  width.  Even  though  we  have 
about  3  laser  photons  per  pixel  to  work  with,  there  are  at  least  189312  pixels  outputs  to 
generate  a  correlation.  Thus,  these  experimental  conditions  are  reasonable.  Variations  in 
aperture  exposure  time  and  separation  give  additional  experimental  latitude. 

Timing  equipment  to  synchronize  the  camera  with  the  laser  pulses  and  the  target 
mount  will  need  to  be  assembled.  We  plan  to  interface  with  the  ISTEF IP9527  Image  Pro¬ 
cessor  and  real-time  disk  array  system  to  digitize  and  store  the  acquired  data. 


Processing 

The  basic  processing  for  missile  typing  is  to  autocorrelate  the  speckle  pattern 
obtained  from  each  laser  pulse  and/or  take  the  Fourier  Transform  of  the  speckle  pattern  to 
get  a  direct  estimate  of  the  autocorrelation  of  the  complex  valued  target  amplitude.  From 
the  autocorrelation  we  can  measure  the  shape  of  the  main  lobe,  which  is  related  to  the 
shape  of  the  object.  From  the  speckle  Fourier  Transform  we  can  examine  the  object  auto¬ 
correlation.  If  the  target  is  not  changing  its  shape  significantly,  e.g.,  it  is  keeping  its  same 
profile  to  the  line  of  sights,  the  correlations  from  the  various  pulses  can  be  summed  to 
improve  the  overall  signal-to-noise  ratio  and  improve  the  shape  discrimination.  When 
there  are  double  exposures  on  the  same  frame  the  autocorrelation  generates  a  secondary 
set  of  peaks.  Measuring  this  spacing  between  the  secondary  peaks  and  knowing  the  expo¬ 
sure  separation  time  will  give  the  speckle  translation  rate. 
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Initial  processing  would  be  done  on  stored  experimental  data  with  available  ISTEF 
computers  or  with  off-site  computers  in  order  to  refine  the  processing  algorithms  using 
actual  acquired  data.  Subsequently,  real-time  processing  with  the  ISTEF  DVP  (Digital 
Video  Processor)  or  other  computers  accessing  the  RTD  (Real-Time  Disk)  array  would  be 
demonstrated. 

Issues 

There  are  a  number  of  technological  issues  to  be  investigated  as  part  of  the  experi¬ 
ments.  They  include  the  following. 

Spatial  and  temporal  coherence  of  the  Long  Pulse  Laser:  The  coherence  characteristics  of 
the  source  laser  will  need  to  be  known  to  understand  its  impact  on  the  speckle  formation. 
It  is  our  understanding  that  coherence  properties  of  the  long  pulse  laser  under  consider¬ 
ation  have  not  been  experimentally  tested.  Partial  coherence  can  reduce  the  contrast  of  the 
speckle  and  its  desired  information  content.  Thus,  measurements  of  the  lasers  output  sig¬ 
nal,  including  its  temporal  and  spatial  coherence,  will  be  carried  out  early  in  the  program 
in  order  to  help  determine  the  exact  nature  of  subsequent  speckle  experiments. 

Depolarization  of  the  scattered  light  by  the  target:  Depolarization  of  the  laser  light  due  to 
scattering  by  the  target  adds  to  the  complexity  of  the  speckle  pattern.  Polarizers  may  be 
required  at  the  detector  to  enhance  discrimination.  Experience  with  the  targets  of  interest 
will  help  resolve  this  issue. 

Atmospheric  effects  via  the  illumination  beam:  Beam  breakup  can  change  the  shape  of  the 
speckle  and  beam  tilt  can  translate  the  speckles.  It  may  be  necessary  to  monitor  the  atmo¬ 
spheric  conditions  to  keep  track  of  the  severity  of  these  effects  and  relate  the  conditions  to 
the  experimental  results. 
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Atmospheric  effects  via  the  scattered  wavefront:  The  stability  of  the  speckle  pattern  is 
affected  by  the  atmospheric  conditions.  For  single  short  pulse  operation  many  of  the  low- 
level  distortion  atmospheric  time  effects  are  not  significant.  However,  if  severe  atmo¬ 
spheric  distortions  are  typically  present  it  may  be  necessary  to  monitor  the  atmospheric 
conditions  to  keep  track  of  the  severity  of  these  effects  and  relate  the  conditions  to  the 
experimental  results. 

Detector  noise  and  photon  noises:  Noise  limits  the  fundamental  sensitivity  of  the  speckle 
measurement  techniques.  We  expect  these  limitations  will  be  quantified  as  the  experiments 
progress  and  the  level  of  the  noise  sources  is  established. 

Target  skin  scattering  characteristics:  The  target  bidirectional  scattering  profile  will  influ¬ 
ence  the  speckle  size  and  shape  for  3-D  objects  such  as  a  cylindrical  missile.  Specular 
reflections  reduce  the  strength  of  the  speckle  formation  and  glints  could  saturate  the  detec¬ 
tor.  Experiments  with  the  targets  of  interest  will  allow  these  effects  to  be  evaluated  and 
incorporated  as  part  of  the  speckle  analysis. 


Feature  shapes  on  the  target:  The  requirements  for  useful  speckle  discrimination  of  three 
dimensional  shapes  and  features  is  not  clear.  Further  analysis  and  experimental  verifica¬ 
tion  would  be  carried  out. 

Micro-motions  of  the  target  during  the  pulse:  Micro- vibrations  of  the  target  driven  by  the 
mount  motion  or  wind  could  affect  the  speckle  techniques.  We  will  assess  minimizing 
vibrations  during  the  design  of  the  mounts  and  how  the  experiments  are  run.  However, 
such  vibrations  might  be  a  benefit  if  they  help  characterize  the  target  via  their  impact  on 
the  speckle. 
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Effect  of  detector  spatial  and  temporal  responses:  Non-uniformity  of  detector  pixel 
responses  could  affect  measurement  accuracies.  Typically,  these  errors  can  be  corrected 
with  post-detection  processing. 


Multiple  targets  with  different  motion:  Multiple  motion  targets  simultaneously  in  the  field 
of  view  complicates  the  analysis  of  the  speckle  pattern.  We  are  not  initially  addressing  this 
effect  since  we  will  have  a  simple  target  set.  Further  analysis  will  need  to  be  carried  out  to 
understand  the  discrimination  of  the  speckle  techniques  under  these  conditions. 


Linearity  if  ICCD  Camera:  Intensity  non-linearity  response  and  saturation  of  the  camera 
intensifier  may  affect  the  processing  of  the  data.  Typically,  such  point  non-linearity  can  be 
corrected  with  post-processing,  but  saturation  should  be  avoided. 
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7  TRACKING  VIA  SPECKLES 


Introduction 

Three  techniques  are  presented  here  for  tracking  a  laser  illuminated  target  from  its 
speckles.  All  the  techniques  depend  on  measuring  the  speckle  motion  to  infer  target 
motion.  One  technique  uses  two  or  more  short  exposures  of  the  speckles  by  a  detector 
array.  Comparison  (i.e.,  correlation)  of  the  exposures  determines  the  movement  of  the 
speckle  pattern.  The  second  technique  uses  a  single  long  exposure  which  produces  speckle 
streaking  from  which  the  speckle  motion  can  be  determined.  The  third  technique  uses  a 
comparison  between  continuous  time-detection  signals  from  many  spatially  distributed 
detectors  to  establish  the  speckle  motion. 

Double  Exposure  Tracking 

In  double  exposure  tracking  we  take  short  exposures  of  the  speckles  generated  by 
coherently  illuminating  a  translating/rotating  target.  Either  the  technique  of  double  pulsing 
the  laser  or  gating  the  detector  can  be  used  to  produce  two  time-separated  exposures.  Each 
exposure  can  be  a  separate  frame,  or  for  simplicity  of  hardware,  both  exposures  can  be 
recorded  (summed)  onto  the  same  frame  of  the  camera.  The  movement  of  the  speckle  pat¬ 
tern  can  be  measured  from  the  autocorrelation  of  the  of  the  summed  pattern  on  the  frame 
or  by  cross-correlation  of  the  separately  exposed  frames.  If  the  two  speckle  patterns  are 
only  translated  versions  of  the  same  pattern,  then  a  local  correlation  peak  indicates  the 
projected  speed  and  direction  (velocity)  of  the  speckles. 

The  speckle  motion  is  caused  by  relative  angular  change  (rotation)  between  the 
target  and  the  source  and  detectors.  In  general,  the  expression  for  speckle  motion  depends 
on  not  only  the  geometry  and  relative  motions  of  the  scenario,  but  also  on  the  wavefront 
generated  by  the  optics  [7. 1-7.4].  For  far-field  operation  in  a  monostatic  system  where  the 
source  and  detectors  are  co-located  the  speckle  velocity  is  given  by 
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where  O  is  the  instantaneous  angular  rotation  rate  of  the  source/detector  as  seen  from  the 
target,  R  is  the  range  to  the  target,  and  a  is  the  angle  between  the  line-of-sight  and  the  ori¬ 
entation  of  the  rotation  axis.  See  Figure  7-1.  Thus,  we  see  from  the  equation  that  knowing 


Source/ 

Detector 


Q 


Figure  7-1.  Rotation  Geometry. 


just  the  speckle  velocity  does  not  allow  us  to  directly  distinguish  between  the  range,  rota¬ 
tion  rate,  and  tilt  angle.  Additional  information  such  as  range  information  and  trajectory 
dynamics  might  allow  us  to  obtain  the  specific  information,  e.g.,  rotation  rate  or  tilt  that 
we  seek. 

The  movement  of  the  speckles  between  exposures  must  be  less  than  the  aperture 
width  for  any  correlation  to  be  achieved.  We  pick  that  the  total  speckle  translation  be  less 
than  half  the  aperture  width  A  as  given  by 

V speckle (7.2) 

where  AT  is  the  time  between  exposures  and  the  speckle  translation  velocity.  If 

we  use  R-  sina  =  100  km,  Q  =  1  revolutions  per  second,  and  A  =  30  cm,  then  AT  needs  to 
be  less  than  120  nsec. 

Decorrelation  of  the  speckle  pattern  occurs  as  the  target  rotates.  The  decorrelation 
angle  depends  on  the  geometry  and  the  target  shape.  For  a  cylinder  the  decorrelation  angle 
is  given  by  [7.5] 
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e  =  (1/ 

where,  L  is  the  diameter  of  the  cylinder.  Thus,  realizing  that  the  angular  speckle  size  is 
given  by  XlL,  we  see  that  the  decorrelation  distance  d  in  term  of  speckle  widths  is  given  by 


d  _  2e 

^speckle 


(7.4) 


For  L  =  0.76  m  and  A,  =  0.5  ^im  the  above  expression  gives  a  value  of  3487.  We 
would  typically  expect  the  number  of  detectors  in  the  array  be  limited  to  less  than  1000 
and  with  at  least  2  detectors  per  speckle  we  see  that  the  decorrelation  distance  is  more  than 
adequate. 


Speckle  Streaking 

With  the  speckle  streaking  technique,  a  single  long  exposure  is  used  to  record  a 
streaked  speckle  pattern.  Taking  the  Fourier  transform  of  this  streaked  speckle  pattern 
forms  a  one  dimensional  ridge  that  is  perpendicular  to  the  speckle  motion.  Thus,  this  ridge 
points  in  the  direction  of  the  projected  spin  axis  and  variations  of  this  direction  from  frame 
to  frame  provides  information  about  the  angular  dynamics  of  the  target.  This  technique  is 
easier  to  use  than  the  double  exposure  technique  since  it  does  not  require  carefully  timed 
exposures.  Long  exposures  for  linear  speckle  translation  may  smear  the  resulting  speckle 
too  much  to  give  velocity  measurements,  but  will  always  give  the  angular  direction.  Angu¬ 
lar  accuracy  to  better  than  0.1  degree  has  been  demonstrated  [7-6]. 


Temporal  Speckle  Motion  Detection 

Temporal  speckle  motion  detection  uses  a  number  of  fast  detectors  that  capture  the 
intensity  of  the  speckles  moving  past  them.  The  detectors  are  spatially  dispersed  in  an 
advantageous  pattern,  e.g.,  along  a  circle.  By  continuously  cross-correlating  the  signals 
from  detectors  the  direction  and  speed  of  the  speckle  motion  can  be  determined.  Although 
a  filled  array  of  detectors  is  not  needed  for  this  technique,  a  large  number  (probably  >100) 
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of  fast  and  sensitive  (i.e.,  expensive)  detectors  as  well  as  a  fast  correlation  processor  would 
be  needed. 
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ABSTRACT 

The  Super-SVA  algorithm  has  been  shown  effective  for  improving  the  resolution  of  microwave 
synthetic-aperture  radar  (SAR)  images.^  In  this  paper  we  review  the  Super-SVA  algorithm  and 
demonstrate  its  application  to  optical  imaging  systems  using  coherent  laser  illumination.  Examples  of 
such  imaging  systems  include  heterodyne-array  imaging,  holographic  laser  radar,  laser  synthetic- 
aperture  radar,  coherent  phase  retrieval  from  pupil-plane  intensities,  and  SCIP  (sheared-beam 
imaging).  We  consider  the  case  of  a  bright  object  on  a  dark  background.  We  show  the  difference  in 
performance  for  different  features  of  the  image  —  prominent  points  versus  extended,  diffuse  areas. 

Keywords:  Super-resolution,  coherent  imaging,  image  processing,  sidelobes,  bandwidth 
extrapolation,  iterative  algorithm 


1  INTRODUCTION 

The  resolution  of  optical  imagery  is  classically  limited  by  diffraction  from  the  collecting  aperture. 
For  many  years  researchers  have  sought  signal-processing  and  mathematical  approaches  to  recovering 
greater  detail  from  collected  imagery.  Recently  a  new  superresolution  technique,  Super-SVA,^  has 
proven  successful  for  obtaining  improved  resolution  from  microwave  synthetic-aperture  radar  (SAR) 
imagery.  In  part,  it  is  effective  there  because,  at  long  wavelengths,  cultural  objects  are  highly 
specular,  making  point-like  scatterers  common  in  microwave  SAR.  In  this  paper  we  explore  the 
efficacy  of  Super-SVA  to  optical  imagery.  This  would  be  applicable  to  coherent  imaging  modalities 
such  as  heterodyne-array  imaging,  holographic  laser  radar^  laser  synthetic-aperture  radar,^  coherent 
phase  retrieval  from  pupil-plane  intensities,^  and  SCIP  (sheared-beam  imaging),^  among  others.  It  can 
also  be  applied  to  non-optical  imaging  modalities  governed  by  a  coherent  transfer  function,  such  as 
magnetic  resonance  imaging  (MRI).  For  these  imaging  modalities  the  impulse  response  is  the  Fourier 
transform  of  the  aperture  function,  rather  than  the  squared  magnitude  of  the  Fourier  transform  of  the 
aperture  function. 

Even  though  the  coherent  optical  imaging  process  may  be  mathematically  equivalent  to  the  SAR 
imaging  process,  the  superresolution  algorithm  may  not  work  as  well  because  of  the  nature  of  the 
optical  imagery.  In  particular,  at  optical  wavelengths,  objects  are  less  specular,  resulting  in  fewer,  if 
any,  point-like  image  features  that  are  important  to  the  success  of  Super-SVA.  In  this  paper  we  show 
examples  of  simulated  coherent  optical  images  and  demonstrate  the  effectiveness  of  Super-SVA  in 
obtaining  additional  resolution  from  them.  We  show  the  difference  between  objects  having  specular 
components  and  those  that  are  purely  diffuse.  We  begin  by  reviewing  Super-SVA,  then  show  results 
for  the  two  types  of  optical  images,  and  discuss  the  results. 
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2SVA  AND  SUPER-SVA 


Before  describing  Super-SVA,  we  briefly  describe  the  spatially  variant  apodization  (SVA) 
algorithm,^  the  engine  that  makes  it  work. 

Sidelobes  of  the  impulse  response  in  an  image  can  be  controlled  by  placing  a  weighting  or 
apodization  function  over  the  aperture.  However,  while  weighing  down  the  edges  of  the  aperture  does 
decrease  sidelobes,  it  does  so  at  the  expense  of  widening  the  main  lobe  of  the  impulse  response, 
thereby  degrading  resolution.  SVA  circumvents  this  trade-off,  yielding  an  image  of  a  point  source 
having  the  fine  resolution  of  an  unweighted  aperture  while  having  no  sidelobes  whatsoever.  It 
accomplishes  this  feat  by  effectively  computing  the  images  from  a  family  of  aperture  weighting 
functions,  and  selecting,  on  a  pixel-by-pixel  basis,  the  image  value  having  the  minimum  magnitude. 
Because  of  the  properties  of  the  cosine-on-a-pedestal  weighting  functions  —  which  includes 
unweighted,  Hamming,  and  Hanning  functions  as  special  cases  —  this  minimization  can  be  done  very 
simply  and  efficiently  from  the  image  computed  from  an  unweighted  aperture.  The  special  property 
that  makes  this  efficiency  possible  is  the  fact  that  the  discrete  Fourier  transform  of  these  weighting 
functions  consists  of  only  three  non-zero  pixels:  the  center  pixel  and  two  adjacent  sidelobe  pixels,  the 
values  of  which  are  proportional  to  the  coefficient  of  the  weighting  function.  Hence  the  images  from 
all  the  other  members  of  the  family  of  weighting  functions  can  be  computed  by  a  simple  three-point 
convolution.  We  can  do  this  in  analytic  form  and  arrive  at  an  expression  for  the  intensity  of  any  given 
pixel  as  a  function  of  the  weighting-function  coefficient.  We  can  also  analytically  solve  for  the  value 
of  the  coefficient  that  minimizes  the  intensity  at  a  given  pixel.  Hence  we  can  efficiently  select,  on  a 
pixel-by-pixel  basis,  the  weighting  function  that  minimizes  the  sidelobes  at  each  pixel.  The  result  is  a 
preservation  of  the  narrowest  mainlobe  from  the  family  of  weighting  functions  but  the  elimination  of 
sidelobes,^  as  illustrated  in  Figure  1  for  an  image  consisting  of  a  single  point  scatterer. 

We  see  from  Figure  1(c)  that  the  result  of  SVA,  in  eliminating  the  sidelobes  while  preserving  the 
fine-resolution  main  lobe,  is  to  produce  an  image  that  has  discontinuous  first  derivatives.  These 
discontinuities  cause  the  Fourier  transform  of  the  SVA  image  to  have  energy  outside  the  area  of  the 
original  aperture,  and  differ  somewhat  from  the  measured  data  within  the  original  aperture. 

These  ingredients  are  all  analogous  to  those  employed  for  the  Gerchberg  super-resolution 
algorithm.^  For  the  Gerchberg  algorithm  we  assume  that  the  object  being  imaged  is  of  finite  extent, 
so  we  truncate  the  image  to  that  finite  extent,  causing  a  discontinuity  that  introduces  higher-frequency 
data  into  its  Fourier  transform.  Then  we  replace  the  Fourier  transform  of  that  image  with  the 
measured  Fourier  data  within  the  area  of  the  aperture,  while  keeping  the  new  information  at  the 
higher  spatial  frequencies.  We  then  recompute  the  image  from  this  extended  Fourier  transform  and 
repeat  the  process. 

At  its  most  basic  level,  the  Super-SVA  algorithm  is  a  modified  version  of  the  Gerchberg 
algorithm,  replacing  the  image-domain  support-constraint  (truncation)  operation  by  the  SVA 
sidelobe-elimination  operation.  As  iterations  progress,  more  and  more  energy  is  added  at  the  higher 
spatial  frequencies  outside  the  aperture,  thereby  extrapolating  the  aperture,  and  achieving  improved 
resolution  in  the  image.  Additional  details  can  be  found  in  Reference  1.  Super-SVA  is  effective 
because  SVA  injects  discontinuities  throughout  the  image  rather  than  just  at  the  edges  of  the  image, 
thereby  creating  higher  spatial-frequency  information  in  a  highly  effective  manner. 
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Figure  1.  Image  of  a  Point,  (a)  From  unweighted  aperture  (dashed  line)  and  from  Hanning- weighted 
aperture  (solid  line),  (b)  after  SVA. 

3  SUPER-SVA  FOR  OPTICAL  IMAGERY 

Super-SVA  works  well  for  microwave  SAR  imagery,  in  which  bright,  point-like  scatterers  are 
common.  The  question  that  concerns  us  here  is  how  well  might  it  work  for  coherent  optical  imagery, 
for  which  specular  returns  are  less  common,  and  sometimes  are  entirely  absent.  We  answer  this 
question  by  computer  simulation  experiments.  We  start  with  an  incoherent  optical  image  of  an  object, 
simulate  realizations  of  speckled,  coherent  imagery  of  the  object,  run  Super-SVA,  and  examine  the 
results.  Since  subtle  details  are  difficult  to  discern  from  a  single  speckled  image,  we  noncoherently 
average  many  images  together  to  reduce  the  speckle  artifacts  and  more  readily  discern  the  difference 
between  the  images  before  and  after  Super-SVA. 

Figure  2  shows  the  results  for  our  first  test  object,  a  photograph  of  a  model  of  a  satellite.  This 
object,  made  of  optically  shiny  materials,  has  substantial  specular  reflections  as  well  as  a  diffuse 
component.  Figure  2(a)  shows  the  object  model.  We  simulated  36  realizations  of  coherent  speckled 
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images,  each  obtained  by  multiplying  the  square  root  of  the  object  by  a  pure-phase  function  having  a 
different  random  phase,  then  low-pass  filtering  using  a  square  aperture  of  length  32  pixels  in  the 
Fourier  domain.  Individual  speckled  images  are  difficult  to  interpret,  so  we  show  in  each  case  the 
noncoherent  average  of  all  36  speckle  realizations.  The  averaging  reduces  (but  does  not  eliminate)  the 
speckle  artifacts.  Figure  2(b)  shows  the  original  diffraction-limited  averaged  image,  computed  from 
an  unweighted  aperture.  Figure  2(c)  shows  the  averaged  image  when  Taylor  weighting  was  applied  to 
the  apertures.  It  reduces  the  sidelobes,  but  at  the  expense  of  degraded  resolution.  Figures  2(d)  shows 
the  result  of  SVA,  which  has  the  fine  resolution  of  the  image  from  the  unweighted  aperture,  but 
without  the  sidelobes.  Figure  2(e)  shows  the  result  of  Super-SVA,  attempting  a  doubling  of  the 
resolution.  For  comparison  we  show,  in  Figures  2(f),  the  averaged  images  from  an  aperture  of  length 
64  pixels,  giving  true  twice-resolution.  Comparing  2(e)  with  2(b)  and  2(f),  we  see  that  the  Super-SVA 
image  does  have  finer  detail  than  the  original  image;  and  that  in  the  specular  point-like  areas  of  the 
image,  the  finer  detail  matches  that  of  the  true  finer-resolution  image.  For  the  diffuse  portions  of  the 
image,  the  improvements  obtained  by  Super-SVA  are  modest,  and  certainly  do  not  achieve  the  factor 
of  two  improvement  seen  in  the  point-like  areas  of  the  image. 

Figure  3  shows  the  same  things  for  our  second  object,  a  CAD-CAM  image  of  a  missile  that  has  no 
specular  returns.  Again,  SVA  greatly  reduces  the  sidelobes,  thereby  removing  the  haze  in  the  image 
around  the  object.  This  is  a  much  more  difficult  case  for  Super-SVA  or  for  any  other  superresolution 
algorithm.  In  this  case,  the  Super-SVA  image  does  look  somewhat  sharper  than  the  original 
diffraction-limited  image,  but  at  the  price  of  increasing  contrast  of  the  residual  speckles  in  the  image. 
In  particular,  the  point  of  the  missile  in  the  upper-left  comer  appears  sharper  than  in  the  original 
image.  The  details  in  the  area  of  the  fins  in  the  lower  right  comer  appear  slightly  sharper  as  well. 
However,  these  effects  are  subtle  enough  that  they  might  be  difficult  to  discern  in  the  image 
reproduced  on  the  page. 

Figure  4  shows  averaged  cuts  through  the  center  region  of  the  32  noncoherently  averaged  images. 
The  cuts  reveal  the  width  of  the  target.  To  get  a  better  feel  for  the  width,  and  to  further  reduce  the 
effect  of  the  speckles,  we  averaged  over  75  parallel  cuts  through  the  center  of  the  image,  each  shifted 
appropriately  to  line  up  the  center  of  the  missile.  The  top  dotted  curve  shows  the  given  image,  with  no 
aperture  weighting.  The  third-from-top  dashed  curve  shows  the  result  of  SVA,  which  is  considerably 
narrower,  a  marked  improvement.  The  fourth-from-top  solid  curve  shows  the  result  of  Super-SVA,  a 
still  greater  improvement.  For  comparison,  the  second-from-top  dash-dot  curve  shows  the  averaged 
cuts  through  the  true  twice-finer-resolution  image  (no  aperture  weighting),  and  the  lowest  dotted 
curve  shows  the  result  of  SVA  on  the  tme  twice-finer  resolution  image.  The  Super-SVA  algorithm 
reveals  an  apparent  object  width  that  is  more  consistent  with  the  tme  twice-finer-resolution  image. 

4  CONCLUSIONS 

We  show  that  SVA  and  Super-SVA,  originally  developed  for  microwave  synthetic-aperture  radar, 
also  work  well  for  coherent  optical  imagery.  We  found  that,  in  the  diffuse  areas  of  an  image,  SVA 
yields  a  considerable  improvement  in  image  quality  over  a  conventional  image,  and  Super-SVA  yields 
a  slight  additional  improvement.  In  the  specular  areas  of  the  image,  SVA  yields  a  considerable 
improvement  in  image  quality  over  a  conventional  image,  and  Super-SVA  yields  a  considerable 
additional  improvement,  resulting  in  a  factor-of-two  improvement  in  resolution.  Further  testing  will 
be  required  to  determine  the  effects  of  noise. 
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On  the  basis  of  these  results,  we  recommend  the  use  of  Super-SVA  on  a  wide  variety  of  imagery 
from  systems  with  coherent  transfer  functions. 
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Figure  2.  SVA  and  Super-SVA  for  Semi-Specular  (Satellite)  Object,  (a)  Object  model,  (b)  diffraction- 
limited  image  from  unweighted  aperture,  (c)  image  from  Taylor-weighted  aperture,  (d)  result  of 
SVA,  (e)  result  of  Super-SVA,  (f)  image  with  twice-fmer  resolution. 


Figure  3.  SVA  and  Super-SVA  for  Diffuse  (Missile)  Object,  (a)  Object  model,  (b)  diffraction-limited 
image  from  unweighted  aperture,  (c)  image  from  Taylor-weighted  aperture,  (d)  result  of  SVA,  (e) 
result  of  Super-SVA,  (f)  image  with  twice-finer  resolution. 
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Figure  4.  Averaged  Cuts  through  the  Missile  Image  (Log  Plots),  (a)  Over  entire  width  of  object;  (b) 
Of  left  edge  of  object.  Curves  (top-to-bottom)  represent  the  given  image  (— •),  the  true  twice-finer 

resolution  image  ( - ),  SVA  on  the  original  image  ( — ),  Super-SVA  on  the  original  image  ( — ),  and 

SVA  on  the  true  twice-finer-resolution  image  (••*•)• 
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Preface 


This  report  documents  the  results  of  a  study  -  Imaging  Correlography  with  Precondi¬ 
tioned  Autocorrelation  Estimates  -  which  was  performed  by  Michigan  Technological  Uni¬ 
versity  for  the  Environmental  Research  Institute  of  Michigan  during  the  period  July  1995 
to  September  1996.  The  principle  investigator  for  this  project  was  Timothy  J.  Schulz  of 
Michigan  Technological  University,  and  the  program  manager  was  Dr.  James  R.  Fienup  of 
the  Environmental  Research  Institute  of  Michigan.  Additional  support  on  the  project  was 
supplied  by  S.  M.  Rathi  of  Michigan  Technological  University. 
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1  Executive  Summary 


Autocorrelation  estimates  obtained  in  imaging  correlography  often  exhibit  large  deviations 
from  the  truth  at  low  lags.  In  particular,  because  the  zeroth-lag  bias  is  calibrated  by 
a  quantity  determined  by  the  data  themselves,  the  zeroth-lag  value  is  often  subject  to  a 
significant  error,  especially  when  few  independent  speckle  frames  are  available.  These  errors 
are  important  because  deviations  in  the  low-lag  autocorrelation  values  from  the  zeroth-lag 
value  are  directly  related  to  a  particular  measure  of  object  smoothness. 

The  contents  of  this  report  include:  (i)  a  brief  review  of  imaging  correlography  and  a 
discussion  of  the  relationship  between  object-estimate  smoothness  and  the  low-lag  values  of 
the  corresponding  autocorrelation  estimate;  (ii)  a  nonlinear  chpping  method  for  precondi¬ 
tioning  the  autocorrelation  estimate  in  an  attempt  to  smooth  (or  regularize)  the  resulting 
object  estimate;  and  (iii)  a  demonstration  of  the  technique  on  laboratory  data. 
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2  Technical  Report 


2,1  Introduction 

Consider  the  imaging  scenario  shown  in  Figure  1.  For  a  fixed  source-object-sensor 
geometry  with  coherent  plane- wave  illumination  of  the  object,  the  complex  amplitude  of 
the  field  in  the  detector  plane,  as  discussed  by  Fienup  in  [1],  can  be  modeled  as 

9t{y)  =  j  (1) 

where; 

•  t  is  an  index  to  time; 

•  X  and  y  are  two-dimensional  spatial  variables  in  the  object-  and  detector- planes,  re¬ 
spectively; 

•  f{x)  is  proportional  to  the  complex  amplitude  of  the  field  scattered  from  the  object 
(in  the  object  plane); 

•  Ot(y)  is  the  time-varying  phase  error  introduced  by  atmospheric  turbulence  (using  the 
single-layer  phase-screen  model); 

•  d  is  the  distance  from  the  object  plane  to  the  detector  plane;  and 

•  A  is  the  nominal  wavelength  of  the  illumination  source. 

This  is  the  Fraunhofer  far-field  propagation  model  [2]  -  quadratic  phase  terms  associated 
with  the  Fresnel  model  have  been  ignored.  For  fully-developed  laser  speckle  [3],  the  far- 


2  B-5 


field  complex  amplitude  gt{y)  is  modeled  as  a  zero-mean  circular  Gaussian  process  with  the 
covariance  function 


E  [9t(y)gKy')]  =  J o{x)e-^^^y-^'y=^dx 

=  ^My)^-jet{y')Q  ,  (2) 

where  o(*)  is  proportional  to  the  incoherent  scattering  function  for  the  object  and  (9(-)  is 
its  two-dimensional  Fourier  transform.  Because  o(-)  is  the  far-held  power-spectral  density 
in  the  absence  of  turbulence-induced  phase  errors,  the  estimation  of  the  object’s  incoher¬ 
ent  scattering  function  o(*)  from  far-held  measurements  is  a  problem  of  two-dimensional 
spectrum  estimation. 


2.2  Object  Estimation  from  Far-Field  Intensity  Measurements 


Measurement  of  the  far-held  complex  amplitude  gt{y)  can  be  accomplished  through 
holographic,  heterodyne,  or  homodyne  methods;  however,  the  short  wavelengths  associated 
with  visible  and  infrared  radiation  make  this  a  difhcult  and  daunting  task.  A  far  simpler 
measurement  process,  and  one  that  is  insensitive  to  turbulence-induced  wavefront  errors,  is 
to  simply  record  the  far-held  intensities: 


m  = 


(3) 


where  jP(-)  is  the  two-dimensional  Fourier  transform  of  the  object-field  complex  amplitude. 
Data  will,  of  course,  be  acquired  only  over  a  limited  spatial  region,  and  if  this  region  is 
defined  by  a  binary  pupil-function  P{y),  then  the  measured  data  can  be  modeled  as 


Dt{y) 


(4) 


If  the  source-object-sensor  geometry  is  slightly  perturbed  over  time,  a  sequence  of  indepen¬ 
dent  speckle  patterns  can  be  obtained.  To  indicate  this,  we  let  Dfc(-)  denote  the  speckle 
intensity  at  time  4,  and  assume  that  the  collection  statistically  independent. 

Image  estimation  from  these  far-field  intensity  measurements  is  called  imaging  correlog- 
raphy  [4-6],  and  most  imaging  techniques  are  motivated  by  the  mathematical  form  for  the 
statistical  moments  of  the  measured  data: 


E[Dk{y)]  =  OiO)Piy), 


(5) 
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and 


=  O^(0)  +  |o(^)|'  P{y)P{y').  (6) 

Because  of  these  relationships,  it  has  been  suggested  in  [4]  that  the  squared-magnitude  of 
the  inverse  Fourier  transform  of  the  intensity  data: 

\dk{yf  =  \j  Dk{y')e=^^^'^' dy'f  ,  (7) 

be  considered  as  the  starting  point  for  an  estimate  of  the  scattering  function’s  autocorre¬ 
lation  function  from  which  a  phase-retrieval  problem  then  be  addressed  for  the  recovery  of 
the  scattering  function  itself.  An  examination  of  the  first  moment  of  this  function  clearly 
reveals  this  motivation: 

E[\dk(y)\^]  =  J I  E[Dk{y')Dk{y")]e^^^y<y’-y"'>dy'dy" 

=  J  J  O\0)P{y')Piy")e^^^y<y'-y'''>dy'dy" 

+ 1  J\o  (^-^)  P{y')P{y'')e^^"^'^^'~^"'>dy'dy" 

=  0^(0)  |p(y)|^  -h  [IpP  +  r^]  (y),  (8) 

where  p(-)  is  the  inverse  Fourier  transform  of  the  pupil  function,  =(=  denotes  the  convolution 
operator,  and 

ro{y)  =  /  1^(1^)  (9) 

is  related  (through  a  scaling  of  the  spatial  variable  y)  to  the  autocorrelation  of  the  ob¬ 
ject’s  incoherent  scattering  function  o(-).  Because  |pp(-)  is  the  incoherent  point-spread 
function  for  a  difFraction-Umited  imaging  system  with  the  pupil  function  P(-),  the  function 
[|pp  *  To]  (•)  is  referred  to  as  the  diffraction-limited  autocorrelation  of  the  object’s  scattering 
function.  Furthermore,  because 

E  y  Dk{y')e^^^^-y'dy'^  f  =  OH0My)\\  (10) 

an  estimate  of  the  diffraction-limited  object  autocorrelation  has  been  suggested  as  [4]: 


J27ry’y\ 


which  has  the  expectation  value: 

.B  [?.(!/))  =  ^[lPp*>-o](s)- 


(12) 
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Although  the  estimate  ro(-)  can  be  used  as  the  data  for  a  phase-retrieval  [7,8]  or  image 
recovery  from  correlations  [9]  problem,  low-lag  values  of  these  data  often  have  a  relatively 
low  signal-to-noise  ratio  (SNR).  As  is  shown  in  the  following  section,  widely  discrepant 
values  at  low-lags  of  an  autocorrelation  estimate  can,  and  usually  do,  lead  to  very  noisy 
(or  rough)  object  estimates.  For  this  reason,  in  the  following  section  we  also  propose  a 
“clipping”  method  whereby  the  autocorrelation  estimate  ro(-)  is  preconditioned  so  that  its 
low-lag  elements  are  close  (or  equal)  in  value. 


2.3  Preconditioning  of  the  Diffraction-Limited  Autocorrelation  Estimate  in 
Imaging  Correlography 


Consider  the  following  measure  of  smoothness  for  a  two-dimensional,  discrete  image: 

5o(A3;)  =  [o(a;)  -  o{x  -  A^,)]^  ,  (13) 

X 

where  Aj,  is  some  predetermined  lag.  If,  for  instance,  Aj,  =  (0, 1),  then  «9o(0, 1)  would 
represent  the  smoothness  as  measured  through  the  intensity  differences  for  adjacent  pixels  in 
the  second  image  dimension.  An  interesting  relationship  between  this  smoothness  measure 
and  the  image’s  autocorrelation  function  is  developed  as  follows: 

<S<,(Aj,)  =  Y  K^)  “ 

X 

=  +  +  J]]o(x)o(x-f  A^,) 

XX  X 

=  2[r,(0)-r„(A,)],  (14) 

where 

‘Toiy)  =  Y  +  y)  (15) 

X 

is  the  autocorrelation  function  for  the  image.  Therefore,  <So(Aa;)  can  be  equivalently  ex¬ 
pressed  through  the  difference  in  values  of  the  image  autocorrelation  at  lags  0  and  Aj,.  If  a 
smoothness  measure  over  various  lags  is  desired,  then  a  global  image  smoothness  might  be 
defined  as 

<5o  =  ^ai<5o(Aa;;),  (16) 

I 

where  a;  is  the  relative  weight  given  to  the  smoothness  in  the  direction  of  the  fth  lag.  This 
quantity  is  then  related  to  the  image  autocorrelation  according  to 


=  2 


L  l 


^o(O)  -  Y 


(17) 


Evidently  then,  images  that  have  small  variations  between  adjacent  pixels  should  also  have 
small  variations  between  low«lag  autocorrelation  values. 

One  method  for  constraining  low-lag  values  of  the  intensity  function  to  be  close  in 
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IT 

Figure  2:  General  pattern  for  the  object  scattering  function  used  to  generate  the  laboratory 
data. 

value  to  the  zero-lag  value  is  to  simply  apply  the  following  nonlinear  processing  to  an 
autocorrelation  function  estimate: 


f  ^o:/ro(j/)/^  o^h  y 

ro{y)=<  I  I  ,  (18) 

I  ro{y),  y^y 

where  3^  is  a  set  of  low-lag  pixel  values  and  {o/}  is  a  set  of  weights  used  for  averaging. 
Another  method  is  to  simply  clip  the  central  values  of  the  autocorrelation  estimate  to  a 
constant  value. 

In  summary, 

the  techniques  described  in  this  section  should  not  be  interpreted  as  methods  for 
improving  the  autocorrelation  estimate^  but  rather  as  methods  for  regularizing 
the  resulting  intensity  estimate.  In  essence,  low-lag  values  of  the  autocorrelation 
estimate  are  used  to  enforce  smoothness  of  the  object  estimate,  whereas  the  large- 
lag  values  are  used  to  recover  the  object  structure. 


2.3.1  Example  with  Laboratory  Data 


Laboratory  data  acquired  by  the  Environmental  Research  Institute  of  Michigan  (ERIM) 
were  processed  by  the  techniques  described  in  the  previous  section.  These  data  were  ac¬ 
quired  by  illuminating  a  pattern  of  retro- reflective  tape  on  glass  in  a  configuration  similar 
to  the  one  shown  in  Figure  2.  The  data  we  processed  represent  data  acquired  over  a  region 
sampled  by  a  128  X  128  uniorm-grid  array.  The  autocorrelation  estimates  obtained  accord¬ 
ing  to  Eq.  (11)  are  shown  in  Figure  3  for  30,  20,  and  10  independent  speckle  realizations. 
One  can  clearly  see  the  problems  with  the  central  lag  in  these  autocorrelation  estimates. 
As  discussed  previously,  this  error  will  lead  to  excessively  noisy  estimates. 


The  autocorrelation  estimates  can  be  modified  by  applying  the  following  thresholding 
operation: 

^max?  ^  ^max 

^0(2/)  ^0(2/)  ^max 


7oiy)  = 


(19) 


where  r^ax  ~  fjm^x{ro{y)}  is  a  percentage  (determined  by  P)  of  the  maximum  value 
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(a)  30  speckle  realizations  (b)  20  speckle  realizations 


(c)  10  speckle  realizations 


Figure  3:  Autocorrelation  estimates  obtained  from  the  laboratory  data  according  to  Eq.  (11) 
for  (a)  30,  (b)  20,  and  (c)  10  speckle  realizations. 

attained  over  the  autocorrelation  estimate  obtained  according  to  Eq.  (11).  The  modified 
estimates  obtained  in  this  way  are  shown  in  Figure  4;  the  scale  factor  used  was:  /?  =  0.10 
for  10  realizations;  /?  =  0.06  for  20  realizations;  and  (3  =  0.07  for  30  reahzations.  The  object 
estimates  obtained  from  these  autocorrelation  estimates  by  using  the  method  proposed  by 
Schulz  and  Snyder  in  [9]  are  shown  in  Fig.  5.  The  object  estimates  have  been  convolved 
with  the  diffraction-limited  point-spread  function  for  the  imaging  system  |pp(-)- 

2.S.2  Issues  for  Further  Consideration: 

Whereas  we  have  established  a  link  between  the  low-lag  autocorrelation  values  and 
object  smoothness,  the  methods  proposed  here  for  preconditioning  the  autocorrelation  esti¬ 
mate  have  not  been  proven  optimal  in  any  sense,  and  other  methods  might  (and  likely  wiU) 
produce  autocorrelation  estimates  from  which  better  object  estimates  can  be  recovered. 
Additionally,  the  use  of  our  clipping  method  requires  the  selection  of  a  clipping  parameter 
/?,  and  a  method  for  selecting  an  optimal  /?  ha5  not  been  addressed.  These  issues  should  be 
considered  in  greater  detail. 


(a)  30  speckle  realizations  (/?  =  0.07)  (b)  20  speckle  realizations  (/?  =  0.06) 


(c)  10  speckle  realizations  (/3  =  0.10) 


Figure  4:  Autocorrelation  estimates  obtained  from  the  laboratory  data  according  to  Eq.  (19) 
for  (a)  30,  (b)  20,  and  (c)  10  speckle  realizations. 


(a)  30  speckle  realizations  (b)  20  speckle  realizations 


(c)  10  speckle  realizations 


Figure  5:  Object  estimates  obtained  from  the  laboratory  data  by  using  the  autocorrelation 
estimates  shown  in  Fig.  4. 


2.3.S  Summary 


In  conclusion,  we  have  presented  in  this  report  a  motivation  and  method  for  precon¬ 
ditioning  autocorrelation  estimates  in  imaging  correlography.  The  main  motivation  lies  in 
the  hypothesis  that  variations  in  the  low-lag  values  of  an  autocorrelation  function  reflect 
the  high-frequency  smoothness  of  the  corresponding  object,  and  if  these  values  are  precon¬ 
ditioned  before  applying  a  phase-retrieval  procedure,  the  resulting  object  estimate  will  be 
smoother  than  one  obtained  from  an  unconditioned  autocorrelation  estimate.  In  addition, 
we  have  presented  a  simple  method  for  preconditioning,  and  demonstrated  the  use  of  this 
method  on  laboratory  data. 
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Appendix  C:  Speckle/Photons  Calculations 

This  MathCad  worksheet  calculates  the  number  of  photons  per  speckle  and  pixel,  and  the  size 
of  the  speckle  as  a  function  of  range. 

Various  parameters  to  be  specified  are: 


X  =0.527'10‘^-m  h  =  6.6256- 10' ^‘’•joulesec  c:  =  3-10' 

E  photon  =  X  ^  photon 

J  =  1 6-joule  Source  Energy  per  Pulse 
Ltx  =0.9-in  Lty  :=ll-m  Target  Dimensions 


,8 


sec 

=  3.7717*  10"’’ • 


•joule 


•t 


Etx'Ety 


D 


xrec 


Dyjec  =  0.0096- m 


^  rec  E)  xrec'^  yrec 


Target  Area 
p  : = 0. 1  Target  Reflectivity 

r  :  =  4-p-A  t  Target  Cross  Section  (Simple  Diffuse  Model)  T  =  3.96 

=  0.0t27-m  Effective  X-Direction  Receiver  Width 
Effective  Y-Direction  Receiver  Width 
Aperture  Area  of  Receiver 
Photon-to-Electron  Conversion  Factor 
Detector  Array  Fill  Factor 
Number  of  Detector  Pixels  in  x  directiion 
Number  of  Detector  Pixels  in  y  directiion 
Atmospheric  One  Way  Transmission  Factor 
Optical  System  Transmitter  Efficiency 
Optical  System  Receiver  Efficiency 
Polarization  Factor 

X-<lirection  Transmit  Beam  Overfill  Factor 
Y-Direction  Transmit  Beam  Overfill  Factor 
Ktx-Ltx 


^  elec  ' 
Kfi„-0.9 

Nx_pix-=768 
Ny_j,ix  =493 
^  atm  = 

’I  trans  ”0-9 
’I  rec  ■  0-9 
hpol  =0.5 
Ktx  =2 

Kty  =2 


Otc(R) 


_  Rty'^ty 


9ty(R)  = 

Sx(R)  =^ 

tx 


X-Direction  Transmit  Beam  Divergence  Angle 

Y-Direction  Transmit  Beam  Divergence  Angle 
X-  Direction  Speckle  size 


A.*R 

S  y( R)  =  j —  Y-  Direction  Speckle  size 

■>'  Lty 

S  x(R)'hl  X  pix 

Q  x(  R)  = - fc - Number  of  detectors  per  speckle  in  X-direction 

^xrec 

S  y(R)'N  y_pix 

Q  y(  R)  = - =: - Number  of  detectors  per  speckle  in  Y-direction 

^  ^  yrec 


Q(R)  -Q  x(R)-Qy(R)  Number  of  detectors  per  speckle 


C-1 


The  photons  per  detector  per  pulse  is  given  by 


Ne(R) 


where  R  is  the  range  to  the  target. 

The  number  of  photons  per  speckle  per  pulse  is  given  by 


N 


spec 


X 

he 


[Ktx-Ltx-(KtyLty)J 


■P’  _  atm  ■’1  trans) 


The  calculated  parameters  for  a  range  of  1  km  are: 

Qy(l -km)  =2.4603  Ne(  1km)  =32.094 


Qx(lkm)  =35.41 
Q(1  km)  =87.1203 


N  spec  =  6.9038- 10^ 


C-2 


speckle  Size  (meters) 


Sx(R) 


0 


X-  Direction  Speckle  size 


XR 


X'R 

S  v(^)  — 0  Direction  Speckle  size 

^  Lty 


Speckle  Size  vs  Range 


The  calculated  parameters  for  a  range  of  1  km  are: 
Sx(  1  km)  =5.8556*  10"“ -m 
Sy(lkm)  =  4.7909*  10“^*m 


D 


xrec 


Sx(l-km) 


=  21.6888 


^yrec 

Sy(lkm) 


=  200.3795 


^  xrec  ^  yrec 
Sx(l-km)  Sy(l  km) 


=  4.346*10^ 
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